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Abstract 

We review the derivation of the metric for a spinning body of any shape 
and composition using hnearized general relativity theory, and also obtain 
the same metric using a transformation argument. The latter derivation 
makes it clear that the linearized metric contains only the Eddington a 
and 7 parameters, so no new parameter is involved in frame-dragging or 
Lense-Thirring effects. We then calculate the precession of an orbiting 
gyroscope in a general weak gravitational field, described by a Newtonian 
potential (the gravito-electric field) and a vector potential (the gravito- 
magnetic field). Next we make a multipole analysis of the potentials and 
the precession equations, giving all of these in terms of the spherical har- 
monics moments of the density distribution. The analysis is not limited 
to an axially symmetric source, although the Earth, which is the main 
application, is very nearly axisymmetric. Finally we analyze the preces- 
sion in regard to the Gravity Probe B (GP-B) experiment, and find that 
the effect of the Earth's quadrupole moment (J2) on the geodetic preces- 
sion is large enough to be measured by GP-B (a previously known result), 
but the effect on the Lense-Thirring precession is somewhat beyond the 
expected GP-B accuracy. 
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1 Introduction 



The Gravity Probe B satellite is scheduled to fly in the year 2000 . It contains 
a set of gyroscopes intended to test the predictions of general relativity (GR) 
that a gyroscope in a low circular polar orbit, with altitude 650 km will precess 
about 6.6 arcsec/year in the orbital plane (geodetic precession) and about 42 
milliarcsec/year perpendicular to the orbital plane (LT precession, see 1^; 
sees. 4.7 and 7.8; [Q, sec. 9.1). In this paper we review the theoretical deriva- 
tion of these effects and in particular consider the contributions of the Earth's 
quadrupole and higher multipole fields. 

We first review the derivation of the metric for a rotating body using the 
standard LGRT approach ( @; g; @, sees. 4.7 and 7.8). The metric is 
characterized by a Newtonian scalar potential (the gravito-electric field) and a 
vector potential (the gravito-magnetic field) (|^, sec. 3.5). We then obtain the 
same result with a simple transformation argument which clarifies the physical 
meaning of the metric (|^, sec. 4.3). Specifically it makes clear that if the 
metric of a point mass contains fundamental parameters such as the Eddington 
parameters a and 7, then to lowest order the metric of a rotating body contains 
no new fundamental parameters Thus there is no new Lense-Thirring or 
frame-dragging parameter to be measured by GP-B or any other experiment. 

We then derive the precession equations for a gyroscope in a general weak 
field system, that is for any scalar and vector potential fields Q . The calculation 
is valid to first order in the fields and velocities of the source body and the 
satellite. The gravitational field of the earth is described by the scalar and 
vector potentials which depend on the shape of the body and the mass and 
velocity distribution inside it. We treat both of these fields by a multipole 
expansion and express the precessions in series of spherical harmonics. We do 
not limit ourselves with the axially symmetric case which has been thoroughly 
studied by Teyssandier (0), ([||). In particular we show that for a solid body 
rotation up to the order I < 2 both precessions depend only on the tensor of 
inertia of the Earth. 

The major nonspherical contributions to the GP-B precessions are from the 
Earth's quadrupole moment. The contribution to the geodetic precession which 
has a magnitude of about 1 part in 10"^ is detectable by GP-B and quite impor- 
tant for the determination of the parameter 7 which is to be measured to about 
1 part in 10^, the most accurate measurement envisioned (Q, sec. 3.5 and in 
particular table 14.2). This contribution has been calculated independently by 
Wilkins Q and Barker & O'Connel ||l^, and then in the most elegant and gen- 
eral way by Breakwell [Tl| ]. The contribution to the Lense-Thirring precession 
is beyond the expected accuracy of GP-B and close enough to the result by 
Teyssandier (0) obtained from a different Earth's model. 

We also estimate the influence of Moon and Sun to find that only the geode- 
tic contribution of the Sun must be included in the GP-B data reduction, as 
anticipated. 
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2 The Metric in Linearized General Relativity 
Theory 

We first briefly review this standard derivation of linearized general relativity 
theory. The metric of a rotating body such as the Earth is obtained by intro- 
ducing a small perturbation of the Lorentz metric r/^j,, that is 

Ofj^L' Vf^i^ + hf^u (1) 

The perturbation /i^^ is assumed to be expressed in isotropic space coordinates 
so that /ill — — /133 — hg. Similarly we suppose the matter producing the 
metric field is described by the energy-momentum tensor of slow-moving and 
low density matter with negligible pressure, 

Tf"" = pu^'u", (2) 

where is the 4-velocity and p is the matter density. The field equations are: 

Gpiy = SnGT^^, (3) 

where G^^ is the Einstein tensor. 

The calculation of G^^ and T^^ to the lowest order in the perturbation is 
straightforward and results in the following equations: 

(^h^^-hr^^,^ ^16nGT^,, h = hl, D^ = -^-V^; (4) 

h\: = (5) 

Here indices are raised and lowered with the Lorentz metric, which is consis- 
tent to lowest order, and the slash denotes differentiation. The last expression 
imposes the so-called Lorentz condition which can always be achieved by a co- 
ordinate transformation (to lowest order) and involves no loss of generality; it 
is also the analog of a gauge choice in electromagnetism. Because we have used 
isotropic coordinates equation (Q) leads immediately to /iqo = hg , and the stan- 
dard classical correspondence gives /loo = hs = 2$. We thereby obtain the wave 
equations for the scalar and vector potentials: 

D^^^-inGp (6) 

D^h^lGirCpv, h = {hoi, ho2, hos}, (7) 

where v is the velocity of the source. These equations may be solved in the 
standard way by means of a retarded Green's function. In the case of a time 
independent system, or a system which changes so slowly that retardation effects 
may be ignored, the solution is 

(8) 



3 



,.«.4G/£<Q|qifi: (9) 



Note that these expressions are analogs of the equations of electrostatics and 
magnetostatics, which is why one may speak about gravito-electric and gravito- 
magnetic effects in the linearized theory. 

In summary, we may write the (Lense-Thirring) line element as 

ds'^ = (l + 2<P)dt^ - (1 - 2$) (if 2 + 2h ■ dfdt (10) 

This is valid to first order in field intensity and source velocity, and will serve 
as a basis for calculating the gyroscope precession. 

Fig. 1 shows the general aspects of the scalar and vector fields of a spinning 
body; the vector field generally points in the same direction as the velocity of 
the surface of the body. 

3 The Metric with Parameters: Derivation by 
Frame Transformation 

It is, in fact, possible to obtain the above result from a different and physi- 
cally interesting perspective, and moreover introduce parameters convenient for 
discussing experimental measurements. Following Eddington, consider the met- 
ric of a massive point with a geometric mass m at a large distance r, so that 
m/r ^ 1. We expand the Schwarzschild solution in isotropic coordinates for 
this situation as 



(1 - m/2r)2 
(1 + m/2r)2 



ds^ = 7^ ' dt^ - (1 + m/2r)^dr'^ 



2m 2m2 \ , 9 / 2m 3m? , 
1 - — + — + . ..jdt^ - n+ — + — + ...]dr' 

Eddington suggested that this be written in terms of parameters as |T^, ||T^ 

, 9 / 2m ^2m2 \ , , 2m \ ,^9 , 

ds^ = [l-a +p +...]dt^ - 1 + 7 + ... df2, (11) 

V r \ r 



where a and (3 and 7 are equal to 1 for general relativity. The power series ( |11[ ) 
is clearly a rather general form for the metric far from a spherical body. 

Since the parameter m which appears in the metric is a constant of integra- 
tion representing the mass of the central body (specifically m = GAf/c^), we 
may absorb the parameter a into it, which is equivalent to taking a = 1. This 
is consistent as long as no independent nongravitational determination of the 
mass of the body is considered. We will nevertheless retain a in our calculations 
as a book-keeping device. 
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Indeed, the parameters in ( |ll| ) may all be viewed as a tool for tracking 
which terms in the metric contribute to a gravitational phenomenon such as 
the perihelion shift of Mercury or the deflection of the starlight by the Sun. 
Alternatively they may be viewed as numbers which may not be equal to 1 if a 
metric theory other than general relativity is actually valid. In either case they 
provide a convenient way to express the results of experimental tests of gravity as 
giving value to the parameters. This parameterized approach has been extended 
to include many other parameters and has been highly developed under the 
name parameterized post-Newtonian theory, or PPN [Q . In this paper we take 
the viewpoint that general relativity is to be tested and emphasize that we are 
not using the more general PPN approach. Solar system measurements give 
(3-1 = (0.2 ± 1.0) X 10-3 and 7 - 1 = (-1.2 ± 1.6) x 10"^, which are of course 
entirely consistent with general relativity. 

We consider now only phenomena in which the quadratic term, ~ m^/r^, 
in goo is unimportant, that is in which we may ignore f3 and assume that the 
underlying gravitational theory is linear. Then for a stationary mass, 

ds^ = ^1 — a ^dt^ — ^1 + 7 ^dr^ (stationary point mass) (12) 

Since this is nearly the Lorentz metric, we may generalize it to a moving mass 
point by simply transforming to a moving system using a transformation that 
is Lorentzian to the first order in velocity: 

tr = t — VX, Xr = X — Vt 

Here the subscript r labels the system in which the mass is at rest, and which 
moves at velocity v in the x direction in the laboratory system. This transfor- 
mation gives the metric for the moving mass point as 



, 1 ( 2m\ , n / 2m\ , ,4m 
ds^ = ( 1 - a jdr - ( 1 + 7 jc^^ + [a + 7) vdxdt, 

(point mass moving in x direction) 
This obviously generalizes for motion in any direction to 



l~a dt^- 1 + 7 Mr 



ds^ = 1 - a dr - 1 + 7 df^ + (a + 7) [v ■ df) dt 

\ r J \ r J r 

(moving point mass) (13) 

As we assume that our theory is linear to this order (like general relativity) , we 
can superpose the fields of a distribution of such point masses and replace m/r 
by $(r) from (^) and 4mv/r by h{r) from (||), resulting in 

p(f') d^r ' 



$(r) = -G j 
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h{f) = 4G 



p{r')v{r')d^r" 



ds^ = (l + 2a$) dt^ - (1 - 27$) dr 2 + (a + 7) (/^ • (14) 

This agrees with the general relativity result (|l2|) when a = 7 = 1 but now 
contains appropriate combinations of Eddington parameters. 

We emphasize that this line element has been obtained for any slowly moving 
mass distribution from the parameterized metric for a stationary mass point by 
transformation and superposition, and thus no new parameter appears in the 
expression. Therefore a measurement of a phenomenon which depends on the 



cross term in (14) provides a value for a + 7 and not for some new parameter 
associated with gravito-magnetism. 

Note also that the result is rather strong since it depends on the observation- 
ally verified Schwarzschild metric ([l^), the well-tested Lorenz transformation, 
and the principle of superposition, valid for any linear theory. 



4 General Precession Relations 

An orbiting gyroscope has its spin axis parallel displaced in accord with the 
metric (Q). We calculate this motion with minimal assumptions about the 
potentials, which could be the potentials of a nearly spherical and rather uniform 
rigid body such as the Earth or a potato-shaped body such as in Fig. 1, with 
some interior mass distribution. We will work always to the first order in the 
potentials and in the velocities, v and V , of the source and orbiting gyroscope. 
The parallel displacement equation for the gyro spin S''^ is ( @; [P, sees. 4.7 
and 7.8; [|], sec. 9.1): 

ds {.^'^ } "'5 

We suppose that the gyro spin 4- vector is perpendicular to the velocity 4- vector, 
which is equivalent to assuming that the gyro spin has no zero component in its 
rest frame. From this assumption the zero component in another frame is easily 
obtained, and to first order in the satellite velocity it is given hy — S -V . 
We calculate the Euler-Lagrange equations in the standard way, and then put 
them into canonical form to give the Christoffel symbols. To lowest order in the 
potential and velocity the Christoffel symbols are 

1 , { i ^ , f « 



ou="*'" loor"''!- i«r-^*i' 
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The Roman indices in these expressions are space indices and run from 1 to 3, 
a slash denotes an ordinary derivative, and I ^ i. Note that the gravitational 

vector potential h occurs in a gauge invariant way, that is only its curl appears. 
Substitution of the Christoffel symbols into the spin equation of motion gives 



5' + (a+7)$|i (S^V^) -jS' -^V' ($1,5') + 



(16) 

We break the drift rate into two parts, the geodetic drift rate due to the scalar 
potential and the Lense-Thirring precession rate due to the vector potential 
h. Separating also symmetric and antisymmetric parts of the geodetic effect, 
we arrive at, in a 3-dimensional vector notation. 



Slt — ^LT X S, 



LT 



a + 7 



V X /i, 



(17) 



and 



SG = nGxs + i-- 



(5-y)v$ + (s- v$)y 

S = Sg + Slti 



+ 7(y • V$)5|, 



(18) 



(19) 



where CIlt and (Ig are the instantaneous values of the Lense-Thirring and 
geodetic precessions, respectively. Since Qlt is the curl of the gravitational 
vector potential the Lense-Thirring precession rate is the analog of the magnetic 
field in magnetostatics theory. The geodetic effect is of the order of the scalar 
potential times the orbital velocity of the satellite, while the Lense-Thirring 
terms are of the order the scalar potential times the velocity of the central 
body (the Earth), which is typically much smaller than the orbital velocity. 
Symmetric geodetic terms are responsible for stretching the spin vector S, and 
the reason for separating them is that their effect almost vanishes when averaged 
over any reasonable satellite orbit. To see this consider the last term in the 
symmetric part of the geodetic drift rate. Using Newton's law in the form 

V$ = —V we may write 




1 dV^ 

2 dt 



2T ' 



where AV^ is the change in the velocity squared in total time T, and it is 
assumed that the drift rate is small and S does not change significantly during 
this period of time. If the orbit is periodic, this quantity will be zero, and for a 
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nearly periodic orbit it will be very small. Similarly for the remaining symmetric 
terms in the geodetic drift we may write: 

( (s . ^) v$ + {s . v$) = - ( {s ■v)^ + is~)v) = 

dt dt / \ dt / T 

Again, this is zero for a periodic orbit and very small for a nearly periodic orbit. 
In summary the average precession rate of the gyro spin is 

0) = 0g) + 0lt), 0g) = (flc) X S, {'Slt) = (fliT) X S, (20) 

with the values of the geodetic and precessions given in ( p7| ) and (p^). 

Note that the geodetic and Lense-Thirring effects are approximately perpen- 
dicular to each other for an approximately polar orbit around a central body 
of reasonable shape; for a circular polar orbit about a spherical body they are 
perpendicular. Fig. 2 shows the orientation of the various vectors for a fairly 
general situation, a satellite in a roughly planar polar orbit about an oddly 
shaped body. The geodetic precession vector is roughly perpendicular to 
the orbit plane and the Lense-Thirring precession vector VIlt has the general 
appearance of a dipole magnetic field. 



5 Effect of Distant Masses 

If there are relatively distant masses, such as the Moon, in the neighborhood of 
the central body, their density distribution poi'T) may be expressed as 



Poir) = ^ MnS{f- fn) 



The effect on the scalar and vector potentials is 



n 

r -rn\ 



It then follows that the precessions due to these distant masses are 



13 ' ^LT = -{(^ + l)2^- 



2 

n . - - . n ' ' " ' 

(21) 

where L„ = M„(r — fvi) x w„ is the angular momentum of the distant mass. 

For the specific case of the Moon and the GP-B experiment the numerical 
values will be discussed in sec. 9. 
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6 Multipole Expansions 



We will henceforth study the case of a body such as the Earth which is rigidly 
rotating, that is v = uj x r. Motivated by expression (|^) we introduce a new 
vector potential quantity n(r): 



h{r) = 4w X n(f), n(r) = G 



p{r')r' dPr' 



(22) 



The vector n(r) is a harmonic function outside the body which we will expand 
in spherical harmonics; for a special spherically symmetric case n(r) is collinear 
with r. In terms of n(r), the metric (Q) may be written as 



ds^ = (l + 2q$) dt^ - (l - 27$) rff ^ + 4(a + 7) (tJ x n) • drdt 



(23) 



Using (H) and (^) it is possible to express the divergence of n(r) via the 
scalar potential, namely, 

v-n = -($ + ?^- V$) 

which allows us to rewrite formula ([T^ ) for the Lense-Thirring precession as 

r^LT = -(0+7)0; 



^LT = -(a + 7) 



(w-V)n+cj $+r-V$ 



n 



z 

(24) 

Here we have chosen the z axis to be along the spin, uj — oj z. 

We now put the origin of spherical coordinates {r, 0, 93} at the center of mass 
of the body, introduce spherical harmonics with the notation 



{C( 
si 



COS mLp, V = c 
sinmiy9, 11 = 3 



and expand the potentials $ and 11 in corresponding series: 



$(f,t) = - 



n,(f,t) 



GM 



l>2. rn, V 



1+ E «LW(-) 



GMR 



.i>i, 



R 



(25) 



(26) 



(27) 



= 1 terms in (^6|) and / = terms in (|27| ) are missing because the origin is 
at the center of mass of the body, see expressions ( |2^ ) and ( p9| ) below]. These 
expansions are in the inertial frame in which the body rotates; the direction of 
the coordinate polar axis z is so far arbitrary; R is the characteristic size of the 
body, which we take to be the equatorial radius for the earth. The potentials 
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are slowly varying functions of time due to the earth rotation, so we write the 
coefficients as explicit functions of time. These coefficients are related to the 
mass distribution in the standard way by 



At) 



_ {2-S„,o) jl-my. 
M {l + m)l 

(2-5,„o) {I -my. 



(28) 



M 



{I 



■),ip)(fr (29) 



Note that generally the domain of integration here is also time dependent. 
In particular, we write for convenience a^g(t) = a;o(i), fl;o(^) = ^ Pwi^) = 
Pl,{t), p-(i) ^ 0. ^ 

The time dependence in ( p8| ) and ( p9| ) may be easily analyzed for a rigidly 
rotating body provided that the rotation axis is fixed both in the inertial space 
and in the body. In that case in the inertial frame with the z axis along the 
spin, uj = Lo z, the earth density is 

P{r, t) = Pe{r, 0,ip- ujt) = pe{r, 6, (p), 

where Peif, 0, (p) = Pe{r) is the time independent density measured in the frame 
rotating with the earth. Substituting this for the density in ( ^ ) and ( ^ 
and transforming to the earth fixed frame we find that the time dependent 
coefficients may be written in terms of constant moments of density according 
to 

(30) 
(31) 

(32) 
(33) 







COS Lut 


— sin Lut 










sin ujt 


cos Lot 






pTmity 




cos Ujt 


— sin ujt 




Plm 


Plm. 




sinujt 


cos Lot 




Plm 



and the time independent coefficients are explicitly given by 
(2 - (5,„o) {I - m) 



Plr. 



M [l + m\ 
(2 - ,5„,o) {I - m)\ 



I, <p) (fr 



M 



{I 



If the rotation axis wanders in the body and/or in the inertial space, the 
relation between the time dependent and time independent coefficients is given 
by a combination of appropriate rotation matrices which is more complicated 
than the one matrix in (|3^) and ( |3l| ) (see c. f. |l^). Generally, a time dependent 
coefficient with the indices I and m is a linear combination of the appropriate 
time independent coefficients with the indices I and n — 0,1, ... ,1. 



10 



Introducing now the time dependent coefficients from ( p2| ) and (|3^) back 
into (|2^) and we obtain the time dependent potentials in a convenient form 
using time independent coefficients only: 



$(f, t) 



n.(r,t) 



GM 



GMR 



l>2. 



.l> 



1, m, z-* 



ifi — Ujt) 



(34) 



(35) 



This is of course what one should expect intuitively; in general this form will 
be most useful for our purposes. 

The constant coefficients af„j in (|^ are those that are measured very ac- 
curately for the Earth (up to / = 18 their values are found in [Q). However, 
for the earth of arbitrary shape and composition it is impossible to express p\'^ 
through a['„j, in other words, their values, and hence the scalar and vector po- 
tentials, are independent. Nevertheless, a useful relationship between the two 
sets of coefficients exists; to describe it, we need a notation for a general moment 
of the density. 



in particular, 



(2-J„,o) {l-m)\ 
M {l + m)\ 



ML 



(36) 



(37) 



Using the definitions ( |32| ) and (33) and the recurrence relations for Legendre 
functions |18 , we derive the following equalities relating pj^^ to a^^^: 

p\!^ = {21 + l)-^[-2-\l + m + 1)(/ + m + 2)ar+i„+i + (2 - <5™i)-^ar+i™-i + 

"(/ - m)\lM{l + m)!] [M^+n-i^+i ~ {I + m - + m)A/r+n_i„,_i] } 
^ = (T)(2/ + + m + 1)(/ + m + 2)a^^,^^^^ + {2- SraiY^ a^+,,^_,~ 

{I - m)\/M{l + m)!] [Ml'^,i_,„,+^ ^ (l + m - 1)(/ + m)M['_^,,_,,^_,] } (38) 

^ = (2? + l)-i{(/ + m + l)ar+i„+i + (2-^,„o)(/-™)!Afr+H-i™A^a + ™-l)!} 

In the second line of (38), the minus sign is taken and fj, — s when v — c, the 
plus sign and fi — c when v = s. 



Pi 
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From ( |34D and (35) using the definition ( |17| ) of J^lt a rather compact mul- 
tipole expansion for the Lense-Thirring precession may be computed: 



a + 7 GMlu 



[{l-m)pl'is^^^-la'j;,J,3 



l>2, m, V 



(39) 

The corresponding expansion for the geodetic precession is too cumbersome to 
be useful. 



7 The Far Field/High Symmetry Approxima- 
tion: Multipoles with / < 2 and the Tensor 
of Inertia 

If the shape of the central body and the mass distribution inside it are known, 
then all the pertinent quantities, including multipole expansion coefficients a"^ 
and may be found by integration, but this is rarely the case. Even when 
af^ are measured, as for the Earth, all the p;^, and the Lense-Thirring effect 
with them, remain entirely undetermined. However, for a body of any shape 
and composition, the coefficients Oj^, {I = 2) andp^"^, {I = 1) can be expressed 
in terms of the elements lij of the tensor of inertia I determined in a standard 
way, 



Writing li, 



020 



li, we find: 

2/3 -h-h 



Pll 



2MR^ 
h 



^22 



h 



-Pio 



Pll 



Pw = 



2A/i?2 



(40) 



Pll 



MR^ ' 



-Pw 



Pll 



'23 



I12 



This is done by comparing the integrals (p: 



(with I = 2) and (H) (with I = 1) 



to lij using explicit expressions of Legendre functions with I = 1,2; formulas for 
020 and are known and used in geodesy for the determination of the Earth's 
moments of inertia |]l6| . 

Introducing (^) into (^), ( |35| ) and dropping the terms with ^ > 2 for $ 
and Z > 1 for n, we first obtain the I < 2 formulas for the potentials: 



$(f) = -- 
r 



I / 3 
M+ — I tri- — (Ir-D 



If-i (tri)? 



(41) 
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From these we then obtain the approximation for the precessions by differenti- 
ation of the expressions (fsl) and (]l7|) (see also (p^)): 



2r3 



M- 



3 / T 5 . ^ . 



{fxV)+—{lrxV) \ (42) 



(tri)- 



(Ir-r) 



uj ■ r 






— (tr l)r— Ir 


UJ + 6 





(43) 



J) and 



[Of course, the same expression for Q.lt is also obtained from 
with 1^2]. 

Note that formulas (40) and ( pl| ) hold in both the body- fixed and inertial 
frames, and the coefficients and moments of inertia in the latter depend on the 
time if the rotation is not axisymmetric, = aj'^(t), = p\'^{t), lij — Iij{t). 
The expressions (^2|) and ( ^^ for the precessions are meaningful in the inertial 
frame where generally I = l{t) = R(t)I(0)R"^(i), R(i) being a rotation matrix 
which converts the body-fixed radius-vector into the inertial one. In the above 
case of a simple rotation about an inertially-fixed axis, uj — luz, 



R(t) 



cosujt —sinujt 
sin cut cos Lot 
1 



The results (|ll)— (||) are valid under either of the two conditions: 

1) far field, R/r < 1; 

2) high symmetry, i.e., all higher order moments are small 

These results alter somewhat our usual notion of the geodetic and Lense- 
Thirring effects: the first one is generally proportional not only to the orbital 
momentum, but to If* x 1^ as well, and the Lense-Thirring precession generally 
points not only in the direction of the angular momentum L = IcD, but has also 
components parallel to uj, f, and Ir. Two particular cases of the inertia tensor 
are of special interest. 

a) Spherical symmetry, I = diag{/, /, /}. In this case, the standard formu- 
las follow immediately from ( ^3|) for a = 7 = 1 : 



3GM 
2r3 



(rxV), 



n 



2GI 



LT 



3 ^: 



(44) 



Note that we have thus shown this to be the exact result for a spherical earth 
with any radial density distribution p — p{r) . 

h) Symmetric top, I — diag{/i, Ii, /}, Ii ^ /. It turns out that for uj = ojz, 
i. e., for the rotation about the material symmetry axis, the previous expression 
for remains true; to lowest order in the oblateness, this also proves to be the 
exact result for a slightly oblate uniform ellipsoid of revolution rotating about 
its semiminor axis. The corresponding expression for Q.g is given and discussed 
in sec. 9. 
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8 Earth Models 



To go beyond the I < 2 approximation, one must make some assumptions about 
the shape of the central body and density distribution inside it, and use the 
data of gravitational potential measurements, as available. Bearing in mind the 
application to GP-B, that is, the Earth, we use the following set of assumptions. 

1. Gravitational potential. We assume it to be measured, i.e., the gravi- 
tational coefficients af^-^ known. Our problem is thus to determine the vector 
potential in terms of the coefficients p]'^^ using as general an earth model as 
possible. 

2. Shape. We assume that the earth is a slightly oblate ellipsoid of revolution, 
so its surface equation to lowest order in the eccentricity e is 

r ^ r,{9,ip) = R{1- ecos^e) (45) 

where R is the semi- major axis (equatorial radius); the eccentricity is the ratio 
of the difference of the semi-axes to the major one. 

3. Mass distribution. We examine two different models: 

a) With pq > and Ap being arbitrary functions of their arguments, we set 

p{r) ^ Po{r) + Ap{e,ip) > 0, J Ap{e, tp) sin dOdip ^ (46) 

sphere 

The first term here describes any depth variation of the average density, and 
the only assumption is that the angular variations are depth-independent. 

b) For arbitrary functions po > and ps, we set 

p{f) ^ Po{r) + Ps{9,ip)S{r ~ rs{0,ip)) > (47) 

Contrary to the previous model, here all angular variations of the density are 
concentrated at the earth's surface. 

We will call model A the set of assumptions 1, 2 and 3o), and model B the 
assumptions 1, 2 and 3b); in both cases the mass distributions are, of course, 
assumed consistent with the gravitational data of the assumption 1. 

As far as the Earth is concerned, assumption 2) refiects the classical Clairaut 
formula (see c. f. |l^) (the eccentricity is e « 3.353 x 10"'^); also the mass model 
36) seems rather realistic for the Earth since the estimated thickness of the layer 
where its mass distribution varies significantly in the angular directions is only 
about 30 km. 

Note that an entirely different Earth model is used in geodesy: it assumes 
that the sum of gravitational and centrifugal potentials is constant at the Earth's 
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ellipsoid's surface (Q)- This allows one to relate aw to the eccentricity and the 
Earth angular velocity only (in particular, to obtain the Earth's gravitational 
oblateness J2 = —020 with a surprisingly good accuracy), but gives zero values 
to a'j^^, m ^ and leaves p]'^ undetermined. 

Evidently, our two models should give bounds for the corrections to the / > 2 
values of the precessions. Moreover, for both of them it is possible to find the 
corrections explicitly by means of the following four steps: 

1° calculate a'j^^ by ( |3^ via spherical harmonics coefficients of the func- 
tion Ap{6, If) for model A, [or Ps{6, for model B]; 

2° since are assumed known, solve the resulting equations for p\^-^^ hence 
having them expressed via aj^^ (thus Ap or ps are uniquely determined at this 
stage through the gravitational data) ; 

3° using that, calculate the needed moments Af^'^-^j -^^jj of the density by 
( ^ in terms of af^; 

4° using the found values of M;';^^j_-^^, express pj^^ through a^^^ according 
to (H). 

In fact, what is described here is a fit of our density distributions ( ^6| ) and 
( p7| ) to the known gravitational coefficients which allows us, in these two 
cases, to express the former through the latter and thus determine uniquely the 
coefficients i.e., the gravito-magnetic part of the field. The implementation 
of this procedure is rather cumbersome, though basically straightforward, so 
the details are given in the Appendix. The results of the calculations to lowest 
order in the Earth's eccentricity are: 



Pirn — 



1 



2{2l + l) 



-{l + m+ + m + 2)a^+i„+i + (? - m){l -m- l)K/a;ii„+i + 



2 — Sml 



^l + lm-1 '^l^l-lm-1 



0{{elf 



Plm 



(T) 



2(2Z + 1) 



Plm 



{l + m+ + m + 2)af_|_i„^i + (/ - m){l - m - l)Kia^_i^+i' 



2 
1 

21 + 1 



2 ( l^ _L A* I 



(48) 



{l + m + l)a;'+i,„ + m)Kia'{_^^ 
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with the minus sign and /i = s for v = c, phis sign and = c for = s in the 
second of these formulas, and 

I + 2 

K; = , model A k; = 1, model B 

£ 4 

It is remarkable that the first non-zero oblateness correction is quadratic, 
and that the results for the two models differ by just a factor (/ + 2/Z + 4) in 
front of a couple of terms. These expressions are needed only for I > 2 since for 
I = 1 formulas ( pO| ) always hold. 



9 Results for Gravity Probe B 

We now apply the above results to the GP-B satellite which will be circling the 
Earth on a low {R/r « 0.9) polar nearly circular orbit. Our aim is to check the 
GR predictions for the geodetic and Lense-Thirring drift of a GP-B gyroscope 
taking into account all the peculiarities of the real Earth's gravito-electric and 
gravito-magnetic field pertinent to the expected experimental error of one part 
in 10^ for the geodetic effect and about one part in 10^ for the Lense-Thirring 
effect. More precisely, one needs to check whether the spherically symmetric 
approximation (^4[), which has been used for many years in the theoretical 
discussion and planning of the experiment, needs any corrections to match the 
above experimental accuracy. 

We use the GR values for the Eddington parameters a = 7 = 1 in the 
sequel. We always use the assumption 2 of the previous section that the Earth 
has a shape of a slightly oblate ellipsoid of revolution with the eccentricity 
e = 3.353 X 10~^. Its semi-minor axis is assumed fixed in the inertial space 
and coincident with the Earth rotation axis and the z axis of the Cartesian 
coordinates. 

1. Gravitational potential. For the required accuracy, it is enough to include 
only the Earth's quadrupolc moment into $, i. e., to set 

- 020 = J2 « 1-083 X 10"^; <,„ = for / = 2, m > 0; I > 2, (49) 

This is because all gravitational coefficients other than the Earth's oblateness 
J2 are at least 2 orders of magnitude smaller ||lj|. By (|49| ) and (^o|), the I < 2 
expression (^0|) for $ is valid with 

I = diag {/i , /i , /}, Ii=I- J2MR^, (50) 

which gives the gravitational potential in the usual form: 

lj>2 

1- J2— P2(C0S(?) 



$(f) 



GM 







GM 










r 
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2. Instantaneous geodetic precession. Just as for the potential, (^9|) assures 
that the I < 2 formula ( p^ ) is valid for D,g- The explicit form of the geodetic 
precession is obtained either from (^2|) and (|5^) or from its general definition 
(Est) by differentiating the above expression for the potential. The result reads: 



3GM 



rxF)+J2?(- 



z 

1-5-2 



fxV]+2z{zxV 



(51) 

This expression contains the J2 corrections to the classical spherically symmetric 
expression (j^). Breakwell obtained it but in a different form, pinning the 
value of the precession to a point on a given satellite orbit rather than to a given 
point in space through which the spinning particle passes with a velocity V. 

3. Gravito-magnetic field and the instantaneous Lense-Thirring precession. 
We now invoke our assumption about the mass distribution (models A and 
B of the previous section) and substitute the values ( ^ ) of the gravitational 
coefficients into the expressions (^^ for p^^. In this way it turns out that for 
both models A and B the only additional to ^ = 1 non-zero coefficients are those 
with ? = 3, m = 0, 1, and their values are 



P31 



3 15 ^ 

^30 = -^ ^2, 



model A; 



P3I 



P3I 



1 7 3 _ 3 ^ 

■ "J 2, P30 — —7; J 2, 



model B; 



(52) 



7 ' 7 

So, for both our models the only addition to the field of a symmetric top is 
the gravito-magnetic field ( |52|) . In terms of the Lense-Thirring precession that 
means that the standard expression (^ ) has to be appended by the terms ( ^2|) 
from the general formula (|39|) for the precession. The result reads: 



GIlu 



LT 



3— - 1 z + 3 ^ 



J2 



GMR^uj R 

.2 1^2 _ ^„„2 n\ 



■Z + A I — ' 



15 ^ 15 ^ , , 
— , An = — , model A; 
98' ° 49' 



where (recall z /r = cos 
Z{cos^ 9) ^ Zo (35 cos"* 61-30 cos^ 9+3), X{coa^ 

3 



(53) 



Zo - — , Xo 



Xo{5cos^9-1); (54) 
3 



7' 



model B 



The second line of ( p3| ) is the correction to the classical Lense-Thirring expres- 
sion in the first line which is induced by the gravito-magnetic field ( |5^ ) . [Recall 
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from section 7, b) , that the dependence of ft lt on J2 for a symmetric top ro- 
tating about its symmetry axis cancels out - remarkably!. Thus the first line of 
( p3| ) contains no J2]. 

Note that the only difference between the Earth's models A and B is a factor 
of 7/5 in the constants Zq and Xq. 

Note also that allowing for more multipoles in the gravitational potential 
than given in the equation ( ^9[ ) we automatically obtain higher order multipoles 
in the gravito-magnetic field by the expressions (48) for p\'^. However, for GP-B 
all these higher order corrections prove to be too small to be taken into account. 



J^. Orbit of the GP-B satellite and the instantaneous values of precessions. 
The perfect orbit for the GP-B satellite would be a circular polar one with the 
altitude hs = 650 km which is described hy r ~ ro = R + hs, V ^ Vq = 
\J GM I ro . In reality it is slightly distorted by the quadrupole moment of the 
Earth's gravitational field; with the lowest order in J2 corrections included, the 
orbit becomes (Q), ([|lO|): 



1 



cos 26* 



(55) 



8 V^o 



- f — 1 sin26'r 
2 V^o 



where f, 0, (p are the corresponding unit vectors. (In view of the orbit symmetry, 
it is enough to consider a half-orbit from one pole to the other, < < tt; 
otherwise, the second half may be assigned to vr < < 27r, with the proper 
direction of the unit vector eg). 

Introducing the orbital velocity ( p5| ) into the formula ( ^ ) and dropping some 
terms 0(J|), we obtain the geodetic precession for the GP-B gyroscope in the 
following form: 



3GAfVb 
2r2 



8 V^o 



. 3 ^ R 
2 V '^0 



1 - 3 cos^ 



Using now the orbit radius (|55| ) and keeping only the corrections linear in J2, 
we arrive at the expression we need: 



BGMVo 



7-0 



5 7 , 

cos 

8 2 



(56) 



To derive a similar formula for JIlt from (g3|) we, of course, need no orbital 
velocity but only the orbit radius (|55|); with accuracy 0(J2) the instantaneous 
value of the Lense-Thirring precession for GP-B turns out to be 
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LT 



o(o) ,0(^2). 



(0) 
LT 



GIlo 



[(3 008^61-1)2 + (3 cos 9 sin 9) x\ 



(57) 



(.-'2) 

LT 



1/ — COS 



9-X 



MR' 



sin 9x 



For convenience, we have assumed that the orbit is in the plane y = 0; the 
values of Z = Z (cos^ 9) and X = X (cos^ 9) are given in (^) , while 



C/(cos^ 



- (6 cos 
4 



5 008^6* + 1), 



y(cos^ 61) = -(2 008^61 - 1) 



(58) 



Note that the expressions ( |56D and ( |57| ) contain the corrections coming 
both directly from the non-spherically symmetric gravito-electric and gravito- 
magnetic field and also through the influence of the former on the orbit. The 
first contribution can be recognized from the second in the Lense-Thirring pre- 
cession ( |5^ ) by the presence of the factor MR^/I, however, in the geodetic 
precession (|6|) this difference is concealed. 

5. Orbit-averaged precessions: the GP-B drift rates. With the above results, 
it is now straightforward to carry out the orbit averaging which reduces to a 
simple integration over 9. In this way, to lowest order in the Earth's oblateness 
J2 « 1.083 X 10-3 using (|^ we find 



r2r 



3GMV 
^rT 



R 



3GMV 



1 - 1.003 X 10" 



(59) 



where the values R = 6,378 km and rg = 7,028fcm for the Earth's and orbit 
radius have been taken. The result exactly coincides with the one found by 
Breakwell [|ll[ in a different way (he gave the formula for the orbit of any 
inclination) . Since GP-B is intended to measure the geodetic precession and the 
Eddington parameter 7 to about a part in 10^, a 0.1% correction is critically 
important. Note that the GP-B gyro spin axis will be initially aligned with the 
reference direction to the Guide Star in the orbit plane, so the geodetic effect 
will cause the in-planc drift of the spin in the North-South direction. 
In a similar fashion, from (57) we derive 



n 



LT 



GIuj 



^ 9 /i?\Vl ry MR^ 



(60) 
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For our model B, Zq = 3/14 according to (| 
for the Earth, 



3); therefore, with MR^/I = 3.024 



LT 



GIu 



8 \ro 



3 Mi?2 
^ " 7 / 



GJg; 



1 - 2.97 X 10" 



(61) 

Thus the J2 correction for the model B is (—0.03%), which is essentially beyond 
the expected GP-B accuracy for the Lense-Thirring effect. For the model A the 
correction is of the opposite sign and almost one order of magnitude smaller, 
about 0.007%. For a different model of the mass distribution inside the Earth 
Teyssandier [Q has obtained a shghtly larger correction of about (—0.011%), 
with the same relevance to GP-B. 



6. Effect of the Moon. According to formula (21 ) from section 5 the geodetic 
effect from a distant mass such as the Moon scales with the mass and inversely 
with the square of distance. Since the Moon has a small mass and is at a large 
distance we expect its effect to be small. A rough estimate for the geodetic 
precession due to the Moon is 



{fla) M \rM 



(62) 



which is too small to be of significance for GP-B experiment. 

The Lense-Thirring effect scales with the angular momentum and inversely 
with the cube of the distance. Since the velocity of the Moon in its orbit is 
small we expect its effect to be small, and the estimate 



Lm 



^ « 10- 



(63) 



shows that it is really far beyond the GP-B accuracy. Thus the moon is of no 
consequence for the GP-B experiment. 

6. Effect of the Sun. Although the Sun is quite distant its mass is large 
compared to the Earth so it is not apparent how large its effect on the precession 
will be. We can calculate this using (^8|). In the inertial frame centered on the 
Sun the velocity of the satellite will be the sum of the satellite velocity in the 
Earth orbit plus the orbital velocity of the Earth 



V^Vs + Veo 



(64) 



But the satellite velocity averages to zero in the course of one orbit, while the 
other factors in (18) change very little, so for long time averages we may neglect 
the satellite velocity compared to the Earth orbital velocity. Treating the Sun 



20 



as a point mass and approximating the Earth orbit as a circle we may thus write 
as 

^ f a + 2-f \ GMsVeo . 

where the unit vector n is perpendicular to the plane of the ecliptic and rs is 
the Earth-Sun distance. The numerical value of this, with a = 7 = 1, is 

(f2g.)Ri 6.3 X — 19 marcsec/yr (66) 

This contribution to the GP-B precession is not negligible and must be included 



in the data analysis. It was first discussed by deSitter in 1916 |21 , p2|, and 
further information and references can be found in the book of Will |4F Note 
an important fact concerning the direction of the above precession vector: it is 
perpendicular to the plane of the ecliptic, whereas the geodetic precession vector 
due to the Earth lies mainly in the equatorial plane of the Earth, so the two are 
not parallel. Indeed the precession due to the Sun is roughly in the direction of 
the Lense-Thirring precession due to the Earth. Finally we may estimate the 
precession due to the spin of the Sun by using (^ ) . Since the spin period of the 
Sun is about 1 month this gives roughly 

(rjf 10~" rad/yr (67) 

which is far too small to be of relevance to the GP-B experiment. 
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Appendix 

Let us show briefly the implementation of the four-step procedure described 
in section 8 which allows us to obtain expressions (^) for gravito-magnetic 
multipole coefficients p]'^. We do it for the Earth's model A; calculations for 
the model B are similar. 

First, we introduce the model A density distribution ( ^6| ) into the definition 
( ^6[ ) of the general moment M^^^ of the density and calculate this moment for 
the shape of a shghtly oblate ellipsoid of revolution (assumption 2 of section 8). 
Working to the first order in the eccentricity e, we find Af^^^^^^ in terms of the 
spherical harmonics coefficients 
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of the function Ap{9,ip) from (p^), namely 



kirn 



k + i 



(Al) 

Here Qim, S";™, Tim arc known positive rational fractions of / and m bounded 
for all pertinent values of those parameters, 

- TO + + TO + 1) {l + m){l-m) 



{2l + l){2l + 3) 
{I -m + l){l -m + 2) 



{2l-l){2l + iy 

{I + m)(l + m - 1) 



(2; + l)(2; + 3) ' (2;-l)(2Z + l) ' 

Formula (A.l) is slightly different for the case I ^ m = 0, but we do not need 
it. For k = I the left-hand side of the equality (A.l) is given via a^,^ according 
to @: 

2 - 6mO {I - my. 

This is a tri-diagonal system of linear algebraic equations for pf^ with small 
off-diagonal elements. Solving it for pf^, we express the latter in terms of a^^: 

'Pin 



1 + 3 
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where 
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l+2m 
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l + l 

{l-m)l 
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7-2m 



0{{df), 



{l + m)\ 

Introducing now the expression (A. 2) back into the general formula (A.l), we 
get all moments, with any k, expressed through the gravitational coefficients: 

^ - ^ ■ " (A3) 



^klm — 



{I + my. 1 + 3 ^ 

■a,, 



_k + 2 
'fc + 3 



2-5mo\{l-myk + i'^^"'' 



with the quantities Qim, Sim, Tim simply related to Qim, Sim, Tim, respectively: 
~ _([+to)! ~ _ {l + m + 2y ^ ~ _(l + rn-2y 



[l - to)! 



{l-m + 2y. 



{I - TO - 2)!^ 



We finally set fc = / ± 1 in (A. 3) and use the resulting expressions in the relation 
(^8|), completing thus the last fourth step of the procedure described in section 
8 and obtaining the equalities (Esl) for pl'^ (model A) . 
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